Antiferromagnetism and superconductivity in layered organic conductors: Variational 

cluster approach 
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The K-(ET)2X layered conductors (where ET stands for BEDT-TTF) are studied within the dimer 
model as a function of the diagonal hopping i' and Hubbard repulsion U . Antiferromagnetism and 
d-wave superconductivity are investigated at zero temperature using variational cluster perturbation 
theory (V-CPT). For large U , Neel antiferromagnetism exists for t' < i^2. with t^2 ~ 0.9. For fixed 
t' , as U is decreased (or pressure increased), a d^2_y2 superconducting phase appears. When U is 
decreased further, the a d^j/ order takes over. There is a critical value of f^i ''^ 0.8 of t' beyond 
which the AF and dSC phases are separated by Mott disordered phase. 
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The proximity of antiferromagnetism (AF) and d-wave 
superconductivity (dSC) is a central and universal fea- 
ture of high-temperature superconductors, and leads nat- 
urally to the hypothesis that the mechanisms behind the 
two phases are intimately related. This proximity is also 
observed in the layered organic conductor k-(ET)2Cu- 
[N(CN)2]C1, an antiferromagnet that transits to a su- 
perconducting phase upon applying pressure [l| (here ET 
stands for BEDT-TTF). Other compounds of the same 
family, k-(ET)2Cu(NCS)2 and K-(ET)2Cu[N(CN)2]Br, 
are superconductors with a critical temperature near 
lOK at ambient pressure. However, another member 
of this family, k-(ET)2Cu2(CN)3, displays no sign of 
AF order, but becomes superconducting upon applying 
pressure[3,l3|. The character of the superconductivity in 
these compounds is still controversial. While many ex- 
periments indicate that the SC gap has nodes (presum- 
ably d-wave), others are interpreted as favoring a node- 
less gap. The literature on the subject is rich, and we 
refer to a recent review article[J| for references. 

The interplay of AF and dSC orders, common to both 
high-Tj, and k-ET materials, cannot be fortuitous and 
must be a robust feature that can be captured in a sim- 
ple model of these strongly correlated systems. k-ET 
compounds consist of orthogonally aligned ET dimers 
that form conducting layers sandwiched between insulat- 
ing polymerized anion layers. The simplest theoretical 
description of these complex compounds is the so-called 
dimer Hubbard model [5|, [g] (Fig. [T]A) in which a single 
bonding orbital is considered on each dimer, occupied by 
one electron on average, with the Hamiltonian 



H = t 
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where Cra (cla) creates an electron (hole) at dimer 
site r on a square lattice with spin projection a, and 
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is the hole number operator, (rr') ([rr']) 
indicates nearest- (next-nearest)-neighbor bonds. As the 
ratio t' /t grows from towards 1, Neel AF is increasingly 
frustrated. The values of t' for the Br, CI, NCS and CN3 
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FIG. 1: (Color online) (A) : Schematic view of the hopping 
terms in the dimer model. V-CPT uses a tiling of clusters 
such as the 4-site cluster drawn here (dashed lines). (B) : 
Qualitative phase diagram inferred from our calculations (not 
to scale). See text for details. 

compounds are thought to be roughly 0.5 — 0.65, 0.75, 
0.75 — 0.85 and 1 respectively [2, |6|]. The local Coulomb 
repulsion U is hard to calculate from first principles, but 
is estimated to be of the same order of magnitude as 
the band-width[J]. We assume that applying pressure 
has the effect of increasing t and t' in proportion, with- 
out affectirig U; thus, for fixed t, we refer to t/U as the 
pressure[5[|6|. Different anions correspond also to differ- 
ent chemical pressures, so that, at ambient pressure, the 
various members of the family of compounds correspond 
to different points on the {U,t') plane. Henceforth we 
will set t to 1, thereby measuring t' and U in units of t. 
The dimer model has been studied with a variety 
of methods. For instance. The fluctuation exchange 
approximation (FLEX) 7] predicts the existence of d- 
wave superconductivity based on diverging susceptibili- 
ties. Quantum Monte Carlo calculations [8] find enhanced 
dSC correlations. Ground state wave- functions obtained 
from variational methods [9j predict a transition between 
a dSC and a spin liquid phase, akin to a Mott transi- 
tion, in the highly frustrated case. The spin liquid phase 
was also predicted by variational methods, and through 
a gauge theory of the Hubbard model [lO|. Variational 



calculations based on the dimer model plus additional 
exchange terms reveal a transition between dSC and an 
AF insulator [ill. Cellular Dynamical Mean Field The- 
ory (COMFxTjial has also been used to reveal the Mott 
transition in that case, without treating the broken sym- 
metry phases. 

We report results of a zero-temperature study of the 
dimer model using Variational Cluster Perturbation The- 
ory (V-CPT)[13]. This method captures short-range 
correlations exactly, treats broken symmetries with a 
rigorous dynamical variational principle, and provides 
dynamical information (the spectral function). It has 
been recently applied with success to the broken symme- 
try phases in models of high- Tc superconductors [l4 1 Iq : 
transitions between AF and dSC states were found as 
a function of doping, both for hole- and electron-doped 
materials, at roughly the correct doping levels. Let us 
summarize our results (see Fig. [TJ3): Beyond the normal 
(disordered) state, the following states have been inves- 
tigated: an AF (Neel) state, a dSC state with (1^2 _y2 
symmetry and an extended s-wave state which, for rea- 
sons explained below, we will call dxy For all values 
of t' we find a (1^2 _y2 phase at low U. For t' smaller 
than a critical value t'^^, we find a dSC to AF transition 
driven by U. For intermediate t' {t'^i < t' < t^j)? the 
AF and dSC phases are separated by a Mott insulator 
phase (MI), with no AF order (spin liquid phase). For 
t' > t'^2: the AF phase disappears, leaving only the dSC 
and MI phases. We find that t'^^ ~ 0.8 and t'^^ ~ 0.9. 
For t' > 0.25 — 0.3, we find a dxy phase at sufficiently low 
U . The possibility of other magnetic orders, in particular 
of a 120° order at t' — t, has not been explored in this 
work, and therefore we make no claim that the MI phase 
corresponds to a spin liquid [3| everywhere. 

Variational Cluster Perturbation Theory. V-CPT is 
an extension of Cluster Perturbation Theory (CPT) 
[ld | based on Potthoff's self-energy-functional approach 
(SFA)[il,[l3. In the SFA one defines a functional f7t[S] 
of the self-energy that is stationary at the physical self- 
energy (here t stands for the matrix tij of one-body 
terms). This functional can be evaluated exactly for a 
self-energy E' that is the exact self-energy of a Hamilto- 
nian H' that differs from the original Hamiltonian only 
through its one-body terms (i.e., H and H' share the 
same interaction part). In V-CPT, the lattice of the 
model is tiled into identical clusters (Fig. [T]) , and H' is 
the restriction of the dimer model to a finite cluster, to 
which one adds various Weiss fields that allow for broken 
symmetries. The Green function of the Hamiltonian H' 
is then calculated exactly (by the Lanczos method) and 
the self-energy functional may then be exactly expressed 
as 
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FIG. 2: (Color online) Pressure it/U) dependence of the an- 
tiferromagnetic (blue triangles), d^2_y2 (red circles, scaled by 
2) and d^j/ (green squares, scaled by 5) order parameter, ob- 
tained with V-CPT on 2 X 2 clusters, for various values of 
i'/i. The vertical lines indicate the transition points, separat- 
ing the various phases. 

where G" is the exact Green function of H' , Gq and Gq 
are the non-interacting Green functions of H and H' re- 
spectively, and fl' is the exact grand potential oi H'. A 
trace over frequencies and wave-vectors of the reduced 
Brillouin zone is carried, and the Green functions carry 
discrete indices related to spin and sites within a clus- 
ter. The functional rJt[E] has become a function rit(t') 
of the parameters of the cluster Hamiltonian H'. The 
task is then to find a stationary point of that function, 
and the exact self-energy of H' at that stationary point, 
denoted S*, is adopted as an approximate self-energy for 
the original model H. 

We insist that the only approximation used is the re- 
striction of the space of self-energies to the set of exact 
self-energies of a family of cluster Hamiltonians, with the 
same interaction part as the original Hamiltonian. In 
particular, short-range correlation effects are treated ex- 
actly through the exact cluster Green function. In order 
to study antiferromagnetism and superconductivity, one 
includes in H' the following "Weiss terms" : 
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n,it') = n'-f ^^indct(i + (G„i-G;,-i)G') (2) 



where Q — {tt,tt) is the Neel wave vector. The dx2-y2 
state is probed by letting Arr' = A if r' = r ± x and 
Arr' — —A if r' = r ± y. These Weiss fields allow for 
the physics of long-range order to seep through the clus- 
ter self-energy. But this is not mean-field theory: the 
interaction term is never factorized, and the Weiss fields 
are not the same as the corresponding order parameters. 
The latter can be calculated from the Green function as- 
sociated with the solution: G*~^ = Gq — T,*. The density 



of electrons n is also computed from the Green function 
G*, as its trace: n = TrG*. 

Gap symmetry. The point group of the dimer model 
is C2v, consisting of a rotation of ir about the z axis, and 
reflections accross the x + y and x — y lines. This group 
admits two possible gap symmetries for singlet supercon- 
ductivity: the dx-^-y-^ state defined above {A2 represen- 
tation), and the isotropic state (s-wave, or Ai represen- 
tation). See Ref. y for a more detailed group-theoretical 
analysis. We also included in our study a SC Weiss fields 
belonging to the Ai representation: 



^r.r+x— y ^r,r— x+y 
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= -A2 



(5) 
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(note that the on-site Coulomb repulsion excludes the 
possibility Ar,r 7^ 0). If Ai = A2, these two possibilities 
form together what is customarily called the dxy state, 
but the relative values of Ai and A2 are not constrained 
by the symmetry of the model. We will nevertheless re- 
fer to this state as dxy Even though this state belongs 
to the Al representation, the corresponding gap function 
will display nodes (roughly along the x and y axes), even 
though these nodes may not be robust - e.g. with re- 
spect to impurities - on group-theoretical grounds [^. 
Recent STM observations [l8j are interpreted as favoring 
this state. 

In this work, five variational parameters were used: the 
Weiss fields M, A, Ai and A2, as well as the chemical po- 
tential fj,' on the cluster, the latter in order to ensure ther- 
modynamic consistency [15|. The cluster chemical poten- 
tial /^', which is part of H', must not be confused with the 
actual chemical potential //, which is part of H and con- 
trols the density of electrons on the lattice. The Nambu 
formalism was used in order to treat anomalous averages. 
Time-reversal violating states such as d -I- is (e.g. a com- 
plex mixture of the Ai and A2 representations) were not 
investigated, because of our use of real-number quantum 
mechanics in numerical computations. 

Results. Fig. [5] shows the order parameters as a func- 
tion of pressure t/U, for different values of t' , using a 
4-site cluster (2 x 2). Fig. [3]shows the same type of data 
obtained on an 8-site cluster (2 x 4). The blue triangles 
were obtained by treating only M and /x' as variational 
parameters, i.e. by forbidding a SC solution. The red 
circles were obtained by treating only A and fi' as varia- 
tional parameters, i.e. by forbidding an AF solution. The 
green squares were obtained by treating Ai, A2 and fi' as 
variational parameters. No mixed solutions were found 
by performing a 3-parameter variation of M, A and fj,' . 
Normal state solutions, where the only variational pa- 
rameter is /i', were also obtained. For each value of U and 
t' , the lattice chemical potential /i was adjusted in order 
for the electron density n to be unity (half-filling) within 
0.1%. Then a comparison of the ground state energy per 
site n + fj,n was made in order to decide which of the 
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FIG. 3: (Color online) Same as Fig. [21 this time for 2 x 4 
clusters. 

four solutions (AF, dx^-y-^, d^y or normal) was preferred 
(only the lowest-energy solutions are shown here). On 
Figs [D and [31 vertical lines separate the different phases 
thus determined, with the SC phases having the lowest 
energy at high pressure, the AF phase at low pressure, 
with an intercalated normal phase - in fact a paramag- 
netic Mott insulator - appearing for t' greater than some 
critical value t'^-^ that is close to 0.8. The AF phase dis- 
appears altogether for t' greater than some critical value 
t'^2 of about 0.9. 

The transitions between the different ordered states 
(AF, dx2-y2 and dxy) are of the first order, whereas the 
transitions from AF to normal look continuous. The re- 
sults obtained on the 4-site and 8-site clusters are in over- 
all agreement on the existence of all four phases and their 
relative location. The critical value Uc{t') for the SC to 
AF/MI transition increases slightly with t' , but stays be- 
tween 5 and 6 for both cluster sizes, at least for t' > 0.75. 
The maximum magnitude of the dx^-y^ order parameter 
does not vary much with t' . By contrast, the critical 
U for the MI to AF transition increases sharply with t', 
as the saturation value of the AF order parameter drops. 
Differences between the two cluster sizes lie mostly in the 
extent of the dxy phase, which is more important in the 
8-site than on the 4-site cluster. This is reflected in the 
"foggy" boundary between the two SC phases that we 
display on the schematic phase diagram of Fig. [1] Note 
that the cluster sizes used are too small to perform any 
kind of finite-size analysis, and we cannot say that the 
results are converged. 

The dxy phase only appears at sufficiently high t' (it 
was not found at t' = 0.25) and low U: it seems favored 
by the dispersion, but at intermediate coupling the tide 
turns towards a dx^-y-^ order. In this V-CPT calculation, 
there is no evidence of a lower critical U: the SC phase 
appears to extend all the way to U — Q and no normal 
metallic phase is found. 

Fig.Hlshows the spectral function for sample solutions 
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FIG. 4: Energy distribution curves (i.e. spectral function) 
for three sample solutions with t' = 0.8t. (A): an antiferro- 
magnetic solution at U = 9f. The SDW dispersion curves 
for the same values of the parameters ([/, fi, t' and t) and 
a mean field h = OAt are also drawn. (B): the same, for a 
Mott insulator solution {U = It). What looked like the AF 
gap in (A) is now a Mott gap. (C): the same, for a d^2_y2 
solution {U — 5t). Notice the d-wave gap maximum along the 
(0, 0) — (0, tt) segment. In all cases, a Lorentzian broadening 
of 0.15i has been used. 

for three different values of U corresponding to (A) the 
AF phase, (B) the MI phase and (C) the dx2-y2 phase at 
t' — 0.8t, on an 8-site cluster. In (A) the mean-field AF 
dispersion is also plotted (full lines), with a mean- field 
h = OAt, a value chosen so as to best match the V-CPT 
dispersions (the actual V-CPT order parameter in that 
case is LM) = 0.497). Note that the mean-field disper- 
sion roughly matches the V-CPT spectral function, at 
least where the weight is important. In (B), the spectral 
function displays a Mott gap across all wave- vectors. It is 
very similar to the AF spectral function (A) , except that 
it does not curve back to follow the AF mean-field dis- 
persion near (0,7r) and (7r/2,7r/2). In (C), we notice the 
dx2-y2 gap, maximum along the (0,0) — (0,7r) segment, 
and vanishing along the (— tt, tt) — (0, 0) segment. Except 
for that dSC gap, the spectral function is that of a metal, 
i.e. there is no other gap in the spectrum. Unfortunately, 



photoemission (ARPES) experiments, which would offer 
a probe on the negative-frequency part of the spectral 
function, have not, to our knowledge, been performed on 
these compounds. 

To conclude, the dimer model of layered organic con- 
ductors captures the essential features of this family of 
compounds: AF and SC phases separated by a first order 
transition, except at higher frustration levels where a MI 
phase (possibly spin liquid) appears. We found that the 
superconducting gap symmetry changes from the Ai (i.e., 
dxy) to the A2 (i.e., d^2_y2) representation of the point 
group as U increases, which leads to the possibility of 
a close interplay between those two states as a function 
of pressure in these compounds. Of course an exciting 
development would be the synthesis of a new member 
of the family displaying Neel order at ambient pressure, 
with a value of t' in the 0.8-0.9 range. Then this com- 
pound would be predicted to go through a spin liquid, 
and then through a dSC phase upon applying pressure. 

We are indebted to A.-M. Tremblay, C. Bourbonnais, 
B. Kyung and D. Fournier for useful discussions. This 
work was supported by NSERC (Canada). Computa- 
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